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ABSTRACT 

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing 
CR propagation models. It was shown that the reacceleration models designed to match secondary to 
primary nuclei ratios (e.g., boron/carbon) produce too few antiprotons. Matching both the secondary 
to primary nuclei ratio and the antiproton flux requires artificial breaks in the diffusion coefficient and 
the primary injection spectrum suggesting the need for other approaches. 

In the present paper we discuss one possibility to overcome these difficulties. Using the measured 
antiproton flux and B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei 
as measured in the heliosphere may contain a fresh local "unprocessed" component at low energies 
perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei 
ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this 
idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni {Z < 28). 
Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances 
of three "radioactive clock" isotopes in CR, ^"Be, ^^Al, '^^Cl, are all consistent and indicate a halo size 
z/i ^ 4 kpc based on the most accurate data taken by the ACE spacecraft. 

Subject headings: diffusion — elementary particles — nuclear reactions, nucleosynthesis, abundances — 
cosmic rays — ISM: general — Galaxy: general 
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1. INTRODUCTION 

The spectrum and origin of antiprotons in CR has been 
a matter of active debate since the first reported detec- 
tions in balloon flights (Golden et al. 1979; Bogomolov 
et al. 1979). Because of the baryonic asymmetry of the 
Universe, antiprotons are not found at rest. There is a 
consensus that most of the CR antiprotons observed near 
the Earth are "secondaries" produced in collisions of en- 
ergetic CR particles with interstellar gas (e.g., Mitchell et 
al. f996). 

The spectrum of secondary antiprotons has a peak at 
about 2 GeV decreasing sharply towards lower energies. 
This unique shape distinguishes antiprotons from other 
cosmic-ray species and allows for searches of primary an- 
tiprotons at low energies. Over the last few years the ac- 
curacy has been improved sufflciently (BESS 1995-2000, 
Orito et al. 2000; Sanuki et al. 2000; Asaoka et al. 2002) 
that we can restrict the spectrum of the secondary compo- 
nent accurately enough to test Galactic CR propagation 
models, and the heliospheric modulation. 



It has been recently shown (Moskalenko et al. 2002) that 
accurate antiproton measurements during the last solar 
minimum 1995-1997 (BESS, Orito et al. 2000) are incon- 
sistent with existing propagation models at the ~ 40% 
level at about 2 GeV while the stated measurement un- 
certainties in this energy range are now ~ 20%. The con- 
ventional models based on local CR measurements, simple 
energy dependence of the diffusion coefficient, and uniform 
CR source spectra throughout the Galaxy fail to reproduce 
simultaneously both the secondary to primary nuclei ratio 
and antiproton flux. 

The reacceleration model designed to match secondary 
to primary nuclei ratios (e.g., boron/carbon) produce too 
few antiprotons because, e.g., matching the B/C ratio at 
all energies requires the diffusion coefficient to be too large. 
The models without reacceleration can reproduce the an- 
tiproton flux, however they fall short of explaining the low- 
energy decrease in the secondary to primary nuclei ratio. 
To be consistent with both, the introduction of breaks in 
the diffusion coefficient and the injection spectrum is re- 
quired, which would suggest new phenomena in particle 
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acceleration and propagation. 

Recently there has appeared some indication that the 
atmospheric contribution to the antiproton flux measured 
in the upper atmosphere is underestimated. If this is true, 
the reacceleration model could still be the best one to de- 
scribe propagation of nuclcon species in the Galaxy (for 
more details see Section 8). However, in this work we have 
assumed the published Galactic antiproton flux, corrected 
for atmospheric production, is accurate. 

In the present paper we discuss another possibility to 
overcome the difficulties encountered by reacceleration 
models. We will show that the inclusion of a local pri- 
mary component at low energies, perhaps associated with 
the Local Bubble (LB), reconciles the data. 

2. INTERSTELLAR COSMIC RAY SPECTRUM 

Just as secondary nuclei, the product of the disinte- 
gration of primary nuclei, are abundant in CR but rare 
in the interstellar medium (ISM), diffuse 7-rays, antipro- 
tons, and positrons are secondary products of interactions 
of mostly CR protons and helium nuclei with interstellar 
gas. The CR propagation model that describes any sec- 
ondary to primary ratio should equally well describe all the 
others: B/C, sub-Fe/Fe, p/p ratios, and spectra of nuclei, 
positrons, and diffuse Galactic continuum 7-rays. 

The diffusive reacceleration models naturally reproduce 
secondary to primary nuclei ratios in CR and agree bet- 
ter with K-capture parent/daughter nuclei ratio (e.g., see 
Jones et al. 2001), though this result is not completely 
conclusive due to the large error bars in CR measurements 
and uncertainties in important isotopic cross sections. It 
is, however, clear that some reacceleration is unavoidable 
in the ISM. Because of the unique shape of the antipro- 
ton spectrum, reacceleration has much weaker effect on it 
than in the case of nuclei. Taking into account that the 
antiproton production spectrum can be calculated accu- 
rately, antiprotons provide a unique complementary tool 
to test propagation models (and heliospheric modulation). 

Our previous result (Moskalenko et al. 2001b, 2002), in 
agreement with calculations of other authors (Molnar & 
Simon 2001), was that matching the secondary /primary 
nuclei ratio B/C using reacceleration models leads to val- 
ues of the spatial diffusion coefficient apparently too large 
to produce the required antiproton flux, when the propa- 
gated nucleon spectra are tuned to match the local proton 
and helium flux measurements. This is an essential short- 
coming. 

Assuming the measured antiproton flux is correct and 

the current heliospheric modulation models are approxi- 
mately right, we have the following possibility to reconcile 
the B/C ratio with the required flux of secondary antipro- 
tons. The spectra of primary nuclei as measured in the 
heliosphere may contain a fresh local "unprocessed" com- 
ponent at low energies thus decreasing the measured sec- 
ondary to primary nuclei ratio. This component would 
have to be local in the sense of being specific to the so- 
^The production cross sections of ^^'^^^N, and 12,13q have been n 



lar neighbourhood, so that the well-known "Local Bubble" 
phenomenon is a natural candidate. 

The idea that CR are accelerated out of SN ejecta- 
enriched matter in superbubbles has been discussed in nu- 
merous papers (e.g., Higdon, Lingenfelter, & Ramaty 1998, 
and references therein). The possibility that the "fresh" 
component is coming from the Local Bubble has been dis- 
cussed by, e.g., MorfiU & Freyberg (1998) and Davis et al. 
(2000). We will hereafter call it the "Local Bubble Hy- 
pothesis." The idea is that primary CR like ^^C and ^^O 
have a local component at low energies, while secondary 
CR like B arc produced Galaxy- wide over the confinement 
time of 10-100 Myr. Then the B/C ratio will be lower at 
low energies than expected in a uniform model, due to the 
enhanced local C (and the reduced Galactic production of 
B). If this idea is correct then the high-energy part of the 
secondary /primary nuclei ratio plus the measured antipro- 
ton flux at maximum, ^ 2 GeV, can be used to restrict 
the value and energy dependence of the diffusion coeffi- 
cient, while the required contribution of the local sources 
can be derived from the measured secondary/primary nu- 
clei ratio at low energies. 

One additional hint for the possible existence of an "un- 
processed" component is the calculated ratio of ^'^C/^^C 
~ 0.14 at 120 MeV/nucleon (modulation potential 500 
MV), which appears to be a factor of two larger than that 
observed when the propagation parameters are tuned to 
the B/C ratio (Moskalenko et al. 2002). The isotope ^^C 
is almost all secondary, as are Be and B isotopes. Since the 
primary source of ^"^C is ^^O, accounting for as much as 
~ 60% of the total, this may indicate an "over-enrichment" 
of the assumed source abundances in oxygen. If so, the 
"over-enrichment" may be true also for primary carbon, 
but tuning to the observed B/C artificially eliminates the 
excess of lithium, beryllium, and boron. We note that 
the ratio of ^^N/^^0 is, however, correct, and the problem 
with overproduction of ^^C may arise from cross section 
errors ^ (see Appendix A). 

3. THE LOCAL BUBBLE HYPOTHESIS 

The low-density region around the sun, filled with hot 
H I gas, is called the Local Bubble (e.g., Sfeir et al. 1999). 
The size of the region is about 200 pc, and it is likely 
that it was produced in a series of supernova (SN) explo- 
sions. Most probably its progenitor was an OB star asso- 
ciation. Though people discuss different scenarios (e.g., 
Maiz-Apellaniz 2001; Berghofer & Breitschwerdt 2002), 
the LB age and the number of SN progenitors appears 
to be similar, ^ 10 Myr and ^ 10 — 20 SN, respectively. 
Most probably they exploded as core-collapse SN II or 
thermonuclear SN Ib/c with a mass of pre-SN stars be- 
tween several and ^ IOA/q, with the last SN explosion 
occuring approximately 1-2 Myr ago, or 3 SN occuring 
within the last 5 Myr. 

There is also some evidence of a SN explosion nearby. 
An excess of ®°Fe measured in a deep ocean core sam- 

:ed only in a narrow energy range. 
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pie of ferromanganese crust suggests the deposition of SN- 
produced iron on earth (Knie et al. 1999). The enhanced 
concentrations were found in two of three layers corre- 
sponding to a time span of < 2.8 Myr and 3.7-5.9 Myr, re- 
spectively. The study suggests a SN explosion about 5 Myr 
ago at 30 pc distance. Another study reports an enhance- 
ment in the CR intensity dated about 40 kyr ago (Sonett, 
Morfill, & Jokipii 1987), which is interpreted as the pas- 
sage across the solar system of the shock wave from a SN 
exploding about 0.1 Myr ago. Taking into account possi- 
ble errors of all these estimates, they point to a nearby SN 
explosion some 1 Myr ago (see also discussion in Benitez, 
Maiz-Apellaniz, & Canelles 2002). 

It could also be that "fresh" LB contributions from con- 
tinuous acceleration in the form of shock waves (Bykov 
& Fleishman 1992), and/or energetic particles coming di- 
rectly from SN remnants still influence the spectra and 
abundances of local CR. The elemental abundances of 
the low-energy nonthermal component in a superbubble 
can differ strongly from the standard cosmic abundances 
(Bykov 2001) due to ejection of matter enriched with 
heavy elements from SN and stellar winds of massive stars 
(Wolf-Rayet, OB stars). The continuous acceleration is 
connected with the lifetime of a shock wave in the LB. A 
reasonable estimate is given by the sound crossing time, 
approximately 2 Myr, for a distance of 200 pc in a lO*' 
K plasma (Bcrghofer & Breitschwerdt 2002). On the 
other hand, the particle crossing time can be estimated 
as t x'^/D ~ 1 Myr for a typical value of the diffusion 
coefficient in the ISM D ^ 10^^ cm s"^ and x ~ 200 pc. 
Therefore, accelerated particles are expected to be present 
in this region. 

4. THE CALCULATION PROCEDURE 

In our calculations wc use the propagation model CAL- 
PROP^ as described elsewhere (Strong & Moskalenko 
1998; Moskalenko et al. 2002); for the present purpose 
the 2D cylindrically symmetrical option is sufficient. For 
a given halo height Zh the diffusion coefficient as a func- 
tion of momentum and the reacceleration or convection pa- 
rameters is determined by data on secondary-to-primary 
ratios in CR. The spatial diffusion coefficient is taken as 
Dxx = PDq{p/ pqY] the corresponding diffusion in momen- 
tum space and other details of the models can be found 
in our earlier papers. We use our standard methodologies 
and include two types of cosmic ray sources. 

Antiproton production and propagation are calculated 
as described in Moskalenko et al. (2002). 

The nucleon injection spectrum of the Galactic compo- 
nent was taken as a modified power-law in rigidity (Jones 
et al. 2001), dq{p)/dp oc p'^ / ^/l + {p/2)-'^, for the in- 
jected particle density. The proton and He spectra are 
tuned to the local measurements as described in Moska- 
lenko et al. (2002). 

The LB spectrum is taken to have the form (as sug- 
gested by Bykov & Fleishman 1992 and Bykov 2001 for 
continous acceleration by interstellar shocks): df /dp oc 

^GALPROP model including software and data sets is available at 



p~^ exp(— p/ pb), where p is the rigidity, and pb is the cut off 
rigidity, ph and the LB source abundances are adjustable 
parameters. In terms of kinetic energy per nucleon E this 
can be re-written as 



m) 



p ''exp(-p/p6). 



(1) 



where a{Z,A) is the abundance of a nucleus {Z,A), Z,A 
are the nucleus charge and atomic number, correspond- 
ingly, m is the atomic mass unit, p is the momentum per 
nucleon, pf, = ^^/{Eb + — m'^- The particular spec- 
tral shape of the LB component is not important as long 
as it decreases sharply towards high energies and is much 
softer than the Galactic CR spectrum. We show the re- 
sults obtained with rj = 1, but they are very similar with 
7] = 2 with pb adjusted correspondingly. 

The procedure to tune the CR elemental abundances, 
secondary /primary nuclei ratios, and antiproton flux we 
adopted was as follows. The high energy part of B/C ra- 
tio and antiproton flux measurements are used to restrict 
the value of the diffusion coefficient and its energy depen- 
dence, while the low energy part of B/C ratio is used to 
fix a value for the reacceleration level and define the pa- 
rameters of the LB component. 

The heliosphcric modulation is treated using the force- 
field approximation (Gleeson & Axford 1968). In our pre- 
vious paper (Moskalenko et al. 2002) we used the best cur- 
rently available model of the heliosphcric modulation, the 
steady-state drift model. The use of the simpler model in 
the current paper is justified by the following arguments. 
First, in the previous paper (Moskalenko et al. 2002) we 
tried to build a model of particle propagation explaining 
both CR nuclei and antiproton fluxes equally well. The he- 
liosphcric modulation calculation was considered as one of 
the main possible reasons for the discrepancy between the 
nuclei and antiproton fluxes. It appears however that the 
interstellar propagation and/or cross section errors are re- 
sponsible for the nuclei/ antiproton discrepancy, unless the 
current models of heliospheric modulation are completely 
wrong, which is not very likely. Second, the current paper 
evaluates a new hypothesis for the origin of CR where the 
main emphasis is on the local CR component. Use of a 
sophisticated modulation model containing additional un- 
known variables will make such evaluation unnecessarily 
complicated. Third, the parameters of the diffusion coef- 
ficient are fixed using the data above few GeV where the 
modulation is weak. 

We consider three different models (Table 1) with pa- 
rameters fixed using the described procedure. They are: 
the simplest plain diffusion model (PD), and two reac- 
celeration models, which use different assumed isotopic 
abundances in the Galactic CR sources and LB sources. 
Diffusive reacceleration model I (DR I) has the same iso- 
topic abundances in Galactic CR and LB sources. Diffu- 
sive reacceleration model II (DR II) is the same as DR I 
except that the LB isotopic abundances are tuned to match 
http://www.gamma.mpe-gaj:ching.mpg.de/~aws/aws.html 
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the low-energy data by ACE and Ulysses, thus increasing 
the freedom to fit the data. The Galactic CR elemental 
source abundances (DR II model) arc tuned to the abun- 
dances measured at high energies where the heliospheric 
modulation is weak. 

The plain diffusion model, without an LB component, 
has already been discussed in Moskalenko et al. (2002). It 
is inconsistent with low energy data on secondary/primary 
ratios, and at high energies matching the B/C ratio would 
cause an overproduction of antiprotons. We do not see 
a plausible modification of this model, even including an 
LB component, which would allow to us simultaneously fit 
antiprotons and the B/C ratio. 

Hence we turn to the models with reacceleration. Fig. 1 
(left) illustrates the process of fixing the normalization of 
the diffusion coefficient using antiprotons. The antipro- 
ton flux shown as calculated in the DR I /II models with 
S = 0.47 and different normalizations in the diffusion co- 
efffcient. Do = 2.6 x lO^^S.S x 10^^,4.3 x lO^* cm s'^ 
at p = 3 GV (for antiprotons this corresponds to kinetic 
energy ^ 2 GeV). The injection index 7 is taken equal 
to 2.28, and the Alfven speed va = 23 km s~^. The an- 
tiproton flux at maximum, 2 GeV, appears to be quite 
sensitive to the value of the diffusion coefficient and allows 
us to fix it at Do = (3.3 ± 0.8) x lO^^ cm s'^. (A la devi- 
ation in the data translates to approximately — 20%-|-25% 
accuracy in Do-) The exact value of S is not critical since 
we compare with the antiproton measurements at max- 
imum, 2 GeV. Inelastically scattered antiprotons, the 
"tertiary" component, appears to be important at low en- 
ergies only in the ISM. Fig. 1 (right) shows corresponding 
calculations of the B/C ratio. A halo height of = 4 kpc 
is used (Strong & Moskakniko 2001; Moskalenko, Mash- 
nik, & Strong 2001a); using values differing by up to a 
factor 2 (the estimated uncertainty) would not affect the 
conchisions, the diffusion coefficient would simply scale ac- 
cordingly. Do oc Zh approximately. In Section 6.2 we re- 
evaluate the radioactive secondaries for the current model 
and show that our adopted value of Zh is in good agree- 
ment with the data. 

The local interstellar (LIS) proton and helium spectra 
used in our antiproton calculations are the best fits to 



the data as described in Moskalenko et al. (2002). The 
proton spectrum is shown in Fig. 2 together with data. 
Because of the measurements with large statistics, mostly 
by BESS (Sanuki et al. 2000) and AMS (Alcaraz et al. 
2000), and weak heliospheric modulation above 10 GeV, 
the error arising from uncertainties in the primary spectra 
is only ~ 5%. The agreement between BESS and AMS 
data, currently considered to be the most accurate, is im- 
pressive. The data collected by other instruments (IMAX, 
CAPRICE, LEAP) at approximately the same solar mod- 
ulation level are lower by 5-10% and have larger error bars. 
Adopting a smaller LIS proton flux would lead to an even 
more dramatic discrepancy between the antiproton flux 
data and the calculation in the "standard" reacceleration 
model (Moskalenko et al. 2002). In any case, allowing for 
a ~ 5 — 10% systematic uncertainty in the proton mea- 
surements does not change the conclusions of the present 
paper because the antiproton data have larger error bars. 

We note that as in the case of other nuclei, there should 
be an LB contribution to proton and He spectra. The 
Galactic injection spectra of protons and He should thus 
be signiflcantly flatter below several GeV to match the 
data points at low energies. This does not influence the 
antiproton production because (i) the LB does not produce 
a significant amount of secondaries, and (ii) the antiproton 
threshold production energy is high, ~ 10 GeV. 

Further tuning can be done using the high energy part 
of the B/C ratio, which is not influenced by heliospheric 
modulation and supposedly contains only a Galactic com- 
ponent of CR. Fig. 3 (left) shows a calculation of the B/C 
ratio for Eb = 500 MeV and different energy dependencies 
of the diffusion coefficient. The plotted curves correspond 
to values of the power-law index d = 0.42, 0.47, 0.52, while 
the injection index was tuned to match the high energy 
spectral data. Index S ^ 0.47 is chosen as giving the best 
match. A pure Kolmogorov spectrum, S = 1/3, thus seems 
to be excluded by antiproton spectrum data taken in com- 
bination with B/C ratio data. 

The B/C ratio as calculated with and without a con- 
tribution of the LB component is shown in Fig. 3 (right). 
The LB component is shown calculated with Ef, = 400, 
500, 600 MeV. It is seen, however, that all three provide 



Table 1 
Propagation parameter sets. 







Diffusion coefficient*^ 


Alfven speed,** 


Source 


Model 


Injection index, 7 


Do, cm^ s 


Index, S 


va/\/w, km s~"* 


abundances 


Plain Diffusion (PD) 


2.16 


3.10 X 10^8 


0.60 






Diffusive 

Reacceleration I (DR I) 


2.28 


3.30 x 10^^ 


0.47 


23 


LBS=CRS 


Diffusive 

Reacceleration II (DR II) 


2.28 


3.30 X 10^8 


0.47 


23 


LBS^^CRS 



Note. — Adopted halo size = 4 kpc. 
^po = 3 GV. 

^va is the Alfven speed, and w is defined as the ratio of MHD wave energy density to magnetic field energy density. 
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Kinetic energy, GeV Kinetic energy, GeV/nucleon 

Fig. 1. — Left: Antiproton flux calculated in DR I/II models with index <5 = 0.47 in the diffusion coefficient and different normalization 
values, Do cm s^^. Solid curves - Do = 3.3 X 10'^'^ at po = 3 GV, upper curve - local interstellar (LIS), lower curve - modulated. Dots 
— Do = 2.6 X 10^**, shown modulated only, dashes — Do = 4.3 X lO"^**, shown modulated only. The two lowest curves (dash-dot) marked 
"tertiary" show separately the LIS spectrum and modulated "tertiary" component for Do = 3.3 X 10^* cm s~^. Modulation was made with 
$ = 550 MV (force field). Data: BESS 95-97 (Orito et al. 2000), BESS 98 (Asaoka et al. 2002) MASS91 (Stochaj et al. 2001), CAPR1CE98 
(Boezio et al. 2001). Right: B/C ratio calculated with LB contribution, Ei, = 500 McV. The lines are coded as on the left panel. Lower 
curves - LIS, upper - modulated (<I> = 450 MV). Data below 200 MeV/nucleon: ACE (Davis et al. 2000), Ulysses (DuVernois, Simpson, & 
Thayer 1996), Voyager (Lukasiak, McDonald, & Webber 1999); high energy data: HEAO-3 (Engelmann et al. 1990), for other references see 
Stephens & Streitmatter (1998). 




10 100 

Kinetic energy, GeV 

Fig. 2. — The proton spectrum as calculated in models DR I/II 
compared with the data (upper curve - LIS, lower - modulated to 
550 MV). Dotted line shows the LIS spectrum best fitted to the data 
above 20 GeV (Moskalenko et al. 2002). Data: IMAX (Menn et al. 
2000), CAPRICE (Boezio et al. 1999), BESS (Sanuki et al. 2000), 
AMS (Alcaraz et al. 2000), LEAP (Sec et al. 1991). 

good agreement with B/C data. Fig. 4 shows the inter- 
stellar and LB carbon spectra in this model. By including 
the LB component, we have therefore been able to obtain 
a model simultaneously fitting p, He, p and B/C data. It 
now remains to apply this model to the full range of CR 
isotopes. 



4.1. Calculation uncertainties 

We do not discuss here possible calculation errors. 

Derivation of such errors is a very complicated matter 
given the many uncertainties in the input, such as the cross 
sections, gas distribution in the Galaxy, systematic errors 
in the CR measurements, heliospheric modulation, atmo- 
spheric corrections etc. Some possible errors and their 
effects have been discussed in Moskalenko et al. (2002). 
Here we qualitatively mention what we think may affect 
our conclusions and what may not. 

Possible errors in the cross section of the abundant CR 
nuclei seem to be less important as they can be compen- 
sated by relatively small corresponding adjustments in the 
primary abundances. They however may be more impor- 
tant in case of less abundant "secondary" nuclei. The cross 
section errors are extensively discussed throughout this pa- 
per. 

Errors in the Galactic gas distribution are not so im- 
portant in the case of stable and long-lived nuclei. Such 
errors are compensated simultaneously for all species by 
the corresponding adjustment of the propagation parame- 
ters (diffusion coefficient). 

Heliospheric modulation may introduce some error, but 
it will be similar for all CR nuclei (or their ratios) because 
the critical parameter here is the charge-to-mass ratio, ap- 
proximately 1/2 for all nuclei except (anti-)protons. 

Systematic measurement errors are difficult to account 
for, but their effect can be reduced by careful choice of the 
data to rely on. This is what we try to do in the present 
paper. 
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Fig. 3. — Left: B/C ratio calculated without (dots) and with (solid, dashes) LB contribution (DR I/II models), E), = 500 McV/nuclcon, 
with different energy dependence in the diffusion coefficient, 5 = 0.52 (short dashes), 0.47 (solid), 0.42 (dashes). Lower curves - interstellar 
(LIS), upper - modulated (<I> = 450 MV). Data as in Fig. 1 (right). Right: B/C ratio calculated without LB contribution (dotted), and with 
Eb = 400 (long dash), 500 (solid), 600 MeV (short dash). Lower curves LIS, upper modulated (^ = 450 MV). Data as in Fig. 1 (right). 

gelmann et al. 1990) at 7.5, 10.6, 16.2 GeV/nucleon com- 
bined with ACE 200 MeV/nucleon data assuming modu- 
lation potential $ = 400 MV. At the chosen HEAO-3 en- 
ergies the heliospheric modulation is weak (for the epoch 
1980 we adopt $ = 800 MV), and it is in the middle of the 
logarithmic interval 0.6 -35 GcV/nucleon covered by the 
HEAO-3 measurements; thus the systematic and statisti- 
cal errors are minimal. To the statistical errors of the ACE 
data, we added 5% systematic error, which is a minimal 
conservative estimate of the uncertainties in the measure- 
ments and modulation potential. 

Fig. 5 shows the quality of the fit. Fig. 5a shows the de- 
viation of the calculated abundances from measurements 
at a given energy expressed in standard deviations. Fig. 5b 
shows the average deviation < S >= ^ ' j- * ? 

where , A™ are the calculated and measured abundances 
for the given energy, and cxj is the standard deviation. 
Fig. 5c shows the relative deviation of the calculated abun- 
dances from measurements at a given energy. 

The DR I fit is systematically low"' at low energies 
(ACE) by as much as 15-30% for elements Z = 11, 13, 19- 
25. It disagrees with high energy abundance of iron by 
~ 20%, and by more than 50% for 27C0, and 28^1. Mean- 
while the high energy data taken separately are consistent 
within 5-10% (Fig. 5c). Because of this low energy dis- 
crepancy we consider further only the model DR II, where 
we allow Galactic CR and LB abundances to be different. 
In this model, the low energy data are used to determine 
the LB source abundances. DR II model provides the best 
fit to all data at the cost of extra free parameters. 

In the DR II model, the Galactic CR source elemental 
abundances are tuned (at a nominal reference energy of 
100 GeV), by a least squares procedure, to the abundances 




Fig. 



1 10 100 1000 

Kinetic energy, GeV/nucleon 

4. — Spectrum of carbon calculated with (solid line, model 



DR II) and without (dashes) LB contribution. Upper curves - LIS, 
lower curves modulated using force field approximation (<P = 450 
MV). Local Bubble (LB) interstellar spectrum is shown by dots. 
Data as in Fig. 10. 



Such effect as the atmospheric correction to the observed 
antiproton flux is very important and may affect our re- 
sults. We discuss it in more details in Sect. 8. 

5. APPLICATION TO NUCLEI UP TO NI 



The DR I model gives an approximate fit to all elements 
(Fig. 5). The source elemental abundances are tuned (at a 
nominal reference energy of 100 GeV), by a least squares 
procedure, to the abundances measured by HEAO-3 (En- 

^This may be due to the errors in the production cross sections employed in the calculations 
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Fig. 6. — Deviation of propagated abundances (model DR II) 
from measured at 4.3, 5.6, 7.5, 10.6, 16.2 GeV/nucleon (HEAO-3, 
Engelmann et al. 1990) given (a) separately for each energy in it's, 
(b) averaged for all five energies in it's, and (c) relative. 



measured by HEAO-3 (Engelmann et al. 1990) at 5.6, 7.5, 
10.6, 16.2 GeV/nucleon (Fig. 6). Fig. 7 shows calculated 
propagated abmidances vs. HEAO-3 data at one particu- 
lar energy 7.5 GeV/nucleon. Relative isotopic abundances 
at the source are taken equal to solar system abundances 
(Anders & Grevesse 1989). The key point in the fitting 
procedure is to obtain the correct abundance of boron. 

Fig. 6 shows the quality of the fit to high energy data, 
where Figs. 6a,b,c show the deviation of the calculated 
abundances from measurements at a particular energy ex- 
pressed in standard deviations, an average deviation, and 
the relative deviation of the calculated abundances from 
measurements at a given energy, respectively. The devia- 
tions from the data at any particular energy are almost all 
within ~ 5%. The calculated abundance of 4Be appears to 
be ~ 7% below the HEAO-3 data. It is, however, the light- 
est nucleus measured by the apparatus, and its measure- 
ments thus may be affected by systematic errors. Some 
disagreement in calculated and measured abundances of 
gF and 23V as seen in Fig. 6b is caused by overproduction 
at only one energy point, 7.5 GeV/nucleon in case of gF 
and 10.6 GeV/nucleon in case of 23V, while at other ener- 
gies calculations agree well with data. 10% overproduction 
of 24Cr is seen only at one energy 10.6 GeV/nucleon. Com- 
pared to the calculations, measurements of ^gK and 27 Co 
are particularly scattered; 27C0 is the least abundant ele- 
ment for Z < 29 and its abundance in CR is measured with 



large error bars. 28-'^i the end of the nucleus charge 
interval measured by HEAO-3 and probably its measure- 
ment is also affected by the systematic errors. In general, 
the deviations in measurements are larger for the least 
abundant nuclei, which is not surprising. 

The LB elemental abundances are tuned simultaneously 
with spectra using the low energy part of the B/C ratio 
and isotopic abundances at 200 MeV/nucleon from ACE 
(Wiedenbcck et al. 2001) and at 185 MeV/nucleon from 
Ulysses (DuVernois & Thayer 1996). For many elements 
ACE and Ulysses abundances differ by 10% (Fig. 8c). For 
this reason, to the statistical errors shown we added 5% 
systematic error. In the same way as for the high energy 
data. Fig. 8 shows the quality of the fit to low energy data 
by ACE and Ulysses. 

The fitting procedure is also influenced by the adopted 
value of the modulation potential. We found that applying 
the following pairs of modulation potentials yields almost 
the same LB elemental abundances: 400 MV for ACE and 
700 MV for Ulysses, 450 MV - ACE and 600 MV - Ulysses, 
and 500 MV - ACE and 500 MV - Ulysses. Other combi- 
nations of modulation potentials make the quality of the 
fit worse. The data are shown to deviate from calculations 
in both directions, which mean that we are unlikely to 
introduce essential systematic error by assuming a wrong 
value of the modulation potential. 
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6. ABUNDANCES IN COSMIC RAYS AND COSMIC RAY 
SOURCES 

6.1. Elemental abundances in cosmic ray sources 

The DR II model with an LB component shows good 

overall agreement with data including secondary to pri- 
mary ratios, spectra, and abundances. The derived Galac- 
tic CR source abundances and LB source abundances are 
given in Table 2 and plotted in Fig. 9 relative to the so- 
lar system abundances. Spectra of boron, carbon, oxy- 



gen, and iron are shown in Fig. 10 for two modulation 
levels, 450 and 800 MV. In case of carbon, the normal- 
ization coefficient in the LB component (eq. [1]) is fixed 
as a(6,12) = 6.35 x 10"'' cm^^ g-i sr^i for r? = 1. The 
calculated sub-Fe/Fe ratio is plotted in Fig. 11. Since the 
elemental abundances are tuned at both high and low en- 
ergies, it agrees well with data. In the intermediate region 
at ~ 1 GeV/nucleon, the agreement could be improved by 
taking Ei, smaller than currently assumed 500 MeV (e.g., 
adopting Eb = 400 MeV may raise it by ~ 10%, similar to 



Table 2 

Elemental abundances? 
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^Normalized to Si=l. 
* Upper limit. 
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the change in the B/C ratio, Fig. 3 right). 

The important result (Fig. 9) is that CR source and LB 
source abundances of all major elements (gC, §0, lo^e, 
i2Mg, uSi, ifiS, isAr, 2oCa, aeFc, asNi) arc in good agree- 
ment with each other. Abundances of the Si group, 2oCa, 
26Fe, 28^1 are near the Solar System abundances. Abun- 
dances of other elements in Galactic CR and LB sources 
are mostly consistent with each other, and with solar sys- 
tem abundances, within a factor of 2. Relative to sili- 
con, gC, 7N, gO, loNe, iqS are underabundant both in 
CR source and LB. This corresponds to the well-known 
first-ionization-potential or volatily correlation (see, e.g., 
Meyer, Drury, & Ellison 1998). Nitrogen in CR sources 
and LB sources differs by a factor of ~ 3, which may be 
connected with production cross section errors affecting 
propagation of Galactic CR component (see discussion in 
Appendix A). 

Secondary nuclei ^gK, 21 Sc, 22^1, 23V appear to be over- 
abundant in the LB sources relative to the solar system 
(shown as upper limits in Fig. 9a), though the derived 
absolute LB abundances are not large. The derived LB 
abundance of 22Ti does not exceed that of 24 Cr while the 
derived abundances of 21 Sc and 23V are not larger than 
that of 25Mn (Fig. 9b). One possible reason for this ex- 
cess is the uncertainty in the production cross sections, 
which is especially large for these nuclei. Sometimes there 



is no measurement at all; in this case one can use only 
phenomenological systematics, which are frequently wrong 
by a factor of two or even more, and/or predictions by 
Monte Carlo codes. Often, there is only one measurement 
at ~ 600 MeV/nucleon, which has to be extrapolated in 
both directions^. This allows only a nearly flat Webber- 
type or Silberberg-Tsao-type extrapolations while the real 
cross sections usually have large resonances below several 
hundred MeV, and decrease with energy above a few GeV 
(see, e.g., Afoskalenko, Mashnik, & Strong 2001a). We 
note that Davis et al. (2000) used semiempirical cross sec- 
tions based on Webber et al. and also predicted fluxes of 
sub-Fe elements which are too low. 

An estimate of the overall error, which is reflected in the 
derived LB source abundances, can be obtained by assum- 
ing the complete absence of ^gK, 2iSc, 22Ti, 23^ 
LB source (shown by crosses in Fig. 9). In this case the 
discrepancy between the calculated propagated CR abun- 
dances of igK, 2iSc, 22Ti, 23^7 ^^'^ thosc measured is be- 
low 20%, and can be removed by allowing the production 
cross sections to increase at low energies by 15 — 20%, 
which seems plausible. 

Another possibility is errors in flux measurements of 
the rare CR species. Fig. 12 shows the calculated abun- 
dances^ tuned at low energies to the ACE and Ulysses 
data. Ulysses and ACE measurements are not always in 



^Only the following reactions are well measured (sec compilation by Mashnik ct al. 1998), on a sample of natural iron consisting mostly of 
®^Fe: p-|-natFe — > 'j^ Sc, 23V, '24 Cr; we use our fits to these data. Reactions producing other isotopes of 2]Sc, 22^'' 23^ 26^'^' ll^'^' 
and 24 Cr have only one or two measurements. There is no data available on the production of 21 Sc, 22'^^! 23^ by 26Fe, 25Mn, and 24Cr isotopes 
(except ||Cr), on production of 2iSc, 22^1. by 23V, and on production of 2iSc by 22'^'- These poorly known cross sections contribute to errors 
on the production of sub-Fe elements at low energies. 

® Calculated Li abundance in the plot shows only secondary lithium produced in CR. 
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Fig. 10. — Spectra of boron, carbon, oxygen, and iron (from top 
to bottom) calculated with LB contribution (model DR II). Upper 
curves - LIS, lower curves - modulated using force field approxima- 
tion ($ = 450 MV - solid curves, * = 800 MV - dashes). Data: 
ACE (Davis et al. 2000, 2001), HEAO-3 (Engelmann et al. 1990), 
for other references see Stephens & Streitmatter (1998) (symbols 
are changed). 



Fig. 11. — Sub-Fe/Fe ratio calculated with LB contribution 
(model DR II), and with E^, = 500 MeV/nuclcon. Lower curve - 
LIS, upper curves — modulated (force field, <E> = 450 and 800 MV). 
Data: ACE (Wiedenbeck et al. 2001; Davis et al. 2000), Ulysses 
(DuVernois & Thayer 1996), Voyager (Lukasiak, McDonald, & Web- 
ber 1997), ISEE-3 (Leske 1993), HEAO-3 (Engelmann et al. 1990), 
Sanriku (Hareyama et al. 1999). 



agrcciincint. Note that even for such an abundant nucleus 
as iron, which is the main contributor to the sub-Fe group, 
the discrepancy exceeds 10%, while the disagreement in 
abundance of 21 Sc is ^ 30%. 

The derived source overabundance of sub-Fe elements in 
the LB could also in principle arise from composition dif- 
ferences between the ISM in the LB and solar or Galactic 
average ISM. This is suggested by the fact that the rela- 
tive abundances of secondary elements in the LB sources 
are systematically larger than in the Galactic CR sources 
(Fig. 9). However, the factors required in case of, e.g., 22Ti 
(Ti/Fe ~ 5% compared to solar or SN 0.1%) appear much 
larger than could reasonably be expected oven for unusual 
SN types. A similar effect may be caused by a specific 
correlation between source and gas density distributions 
(Ptuskin & Soutoul 1990). 

6.2. Isotopic distributions in cosmic rays 

Be and B isotopes arc assumed all secondary, thus there 
is no possibility to tune them. The DR II model calcu- 
lation shows perfect agreement with the data on relative 
isotopic abundances of Be and B (Fig. 13). This is in con- 
trast with a "standard" reacceleration model, where we 
obtained a 15% discrepancy with relative abundances of 
■^Be and ^Be isotopes (Strong & Moskak-nko 2001). 

Abundances of stable isotopes of other elements are not 
very conclusive because they are present in the sources, 



but O and Si isotopic distributions still agree very well 
with data (Webber et al. 1996; Webber, Lukasiak, & Mc- 
Donald 1997; DuVernois et al. 1996; Hesse et al. 1996; 
Wiedenbeck et al. 2001) assuming only ^''O and ^*^Si iso- 
topes are present in the LB component. C and N isotopic 
distributions do not agree too well (Fig. 13), but this may 
point to a problem with cross sections. The calculated ra- 
tio 13C/12C ~ 0.11 at 120 MeV/nucleon ($ = 500 MV) 
in the model with LB contribution is still a factor ~ 1.5 
too large compared to the measured value, 0.0629 ±0.0033 
(Voyager 50-130 MeV/nucleon, Webber et al. 1996) and 
0.078 ± 0.011 (Ulysses 100-200 MeV/nucleon, DuVernois 
ct al. 1996), which may be connected in part with over- 
production of ^^C by ^^N. (A discussion of the cross- 
section uncertainties for C and N isotopes is given in Ap- 
pendix A.) If we replace the cumulative cross section 
^^N-hp — »^^C with cross section ^'^N+p ^^^C, the calcu- 
lated ratio ^^C/^^C will be lowered by 10% as estimated 
(sec Appendix A and Fig. A15). Assuming the absence of 
the isotope ^'^C in the Galactic CR sources gives another 
10% reduction. Altogether these corrections yield ^^C/^^C 
~ 0.09, close to the data. 

Fig. 14 shows calculated ^°Be/^Be, 26A1/27A1, 36ci/Cl, 
^^Mn/Mn ratios, usually used as "radioactive clocks" in 
CR, for a halo size z/j = 4 kpc vs. data. In case of Be, Al, 
CI, the agreement with the most accurate low energy data 
by ACE is very good and all the ratios are consistent with 
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et al. (1998b), besides, it seems too low compared to sim- 
ilar neutron knock out reactions on isotopes of Fe and Cr. 
Reactions ^^Fc+p ^''''^Mn, ''''^Fe+p ^"''"'Mn are not mea- 
sured at all. Similarly, only one ^^Mn major production 
cross section ^^Fe+p — >^^Mn has even one point measured. 



each other indicating that Zh = 4: kpc is a good estimate. 
Higher energy data by ISOMAX (Be) are also consistent 
with calculations considering the large error bars. 

The ^''Mn/Mn ratio indicates a somewhat smaller halo, 
but this may be related to uncertainty in its half-life 
and/or production cross section. The half-life of ^''Mn 
against (3~ decay is the most uncertain among the four 
radioactive isotopes - it is only one which is not mea- 
sured directly. It is derived indirectly based on /3+ de- 
cay branch half-hfe, which yields an estimate ti/2(/3~) = 
6.3± 1.3[stat] ± l.l[theor] x 10^ yr (Wuosmaa et al. 1998). 
In case only the half-life is wrong, to get the ^^Mn/Mn 
ratio consistent with lighter element ratios and with ACE 
data (for Zh = 4: kpc) requires ti/2(/3") ~ 2 Myr (Fig. 14). 
Another estimate of "''"'Mn partial half-life based on CR 
propagation calculations gives ^1 — 2 Myr (DuVernois 
1997). Apart from the half- life, a possible source of er- 
rors can be production cross sections of Mn isotopes. The 
fact that the propagated isotopic abundance of ^^Mn, 
^^Mn/Mn= 0.50, is correct (^^Mn is a K-capture isotope 
that is absent in the ISM) indicates that some important 
production channels of ^'''^^Mn may be not calculated cor- 
rectly. (For instance, in case of stable ^^Mn, we have a 
freedom to c;hoose its abundance in the LB component at 
low energies, which may compensate for underproduction 
of this isotope in CR. Meanwhile, this LB ^^Mn does not 
produce any ^''Mn.) Only the reaction "'''Fe+p ^'"''^Mn 
on a natural sample of Fe has been measured well enough 
(see compilation by Mashnik et al. 1998). The cross section 



7. DISCUSSION 

The proton spectrum at low energies still remains uncer- 
tain. The only secondaries produced below the antiproton 
production threshold are positrons and 7-rays. However, 
the positron spectrum alone cannot provide conclusive in- 
formation on the proton spectrum on a large scale be- 
cause (i) the large energy losses of positrons mean that 
the positron spectrum by its nature is local (ii) and there 
are possibly sources of primary positrons such as pulsars. 
Diffuse 7-rays can provide a tool to test the spectrum 
of protons in distant regions, but the a priori unknown 
contribution of electrons via inverse Compton and brems- 
strahlung complicates the picture. A test of the He spec- 
trum at energies below ~ 10 GeV/nucleon can be made 
using the CR deuteron and "'He measurements similar to 
what was done in this paper for heavier nuclei. We plan 
to address this issue in future work. 

We should mention that there is another possibility to 
get the correct antiproton flux in reacceleration models, 
which is to introduce an additional proton component at 
energies up to approximately 20 GeV. The latter energy is 
above the antiproton production threshold and effectively 
produces antiprotons at ~ 2 GeV and below. The intensity 
and spectral shape of this component could be derived by 
combining restrictions from antiprotons and diffuse 7-rays. 
Interestingly this kind of spectrum was used in our HEMN 
model (hard electrons and modified nucleons. Strong, Mos- 
kalenko, & Reimer 2000) to match the spectrum of diffuse 
7-rays as observed by EGRET (Hunter et al. 1997). The 
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advantage of this approach is that the diffuse 7-rays which 
we observe carry information on the large-scale Galactic 
spectrum of CR (producing antiprotons) while particles 
we measure may reflect only the local region. 

One more (non-standard) interpretation is that the so- 
lar modulation is weaker than assumed, and this would 
eliminate the need for a LB component. (A cornerstone 
of the current theories of hcliosphcric modulation is the 
local interstellar spectrum, which is not known but taken 
a priori.) With a modulation potential as small as ~ 200 
MV one can obtain a consistent reaccclcration model com- 
bining B/C, antiprotons, and other species simultaneously. 
To get an agreement with nucleon spectral data, the in- 
jection spectra in such a model should be flatter at low 
energies than the usually adopted power-law in rigidity. 

Recently there has appeared some indication that the 
atmospheric contribution to the antiproton flux measured 
in the upper atmosphere is underestimated. Monte Carlo 
simulations of the hadron cascade development in the up- 



per atmosphere have shown that the antiproton flux in- 
duced by ]3j4-rcactions on air nuclei is larger, at least, by 
~30% (Huang, Derome, & Buenerd 2001) compared to of- 
ten employed calculations with analytical production cross 
sections. This means that the flux of antiprotons in CR 
in reality may be lower at the top of the atmosphere by 
at least 25-30%. If the latter is true, the reacceleration 
model (even without LB) could still be the best one to de- 
scribe propagation of nucleon species in the Galaxy. The 
inclusion of all known eff'ects such as, e.g., sub-threshold 
antiproton production on the abundant atmospheric N and 
O, may be important for evaluation of the correct atmo- 
spheric background. 

We note that Donato ct al. (2001) claim to have ob- 
tained agreement with antiproton measurements in a reac- 
celeration plus convection model using the parameters de- 
rived from B/C and sub-Fc/Fc ratios (Maurin ct al. 2001). 
Apart from having one more free parameter (convection 
plus reacceleration), they in fact fitted B/C and sub-Fe/Fe 
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ratios only at high energies (since in their fitting procedure 
the high energy data overweights few low energy points). 
Their calculated ratios at low energies arc higher than the 
Voyager and ACE data by approximately 20% or about 6a 
(see their Figs. 3, 4 in Maurin et al. 2001). This is, how- 
ever, where most of the problem lies. Besides, in their nu- 
clear reaction treatment, they use the semiempirical cross 
sections by Webber et al., which are not particulary accu- 
rate at low energies while rcaccclcration models are sensi- 
tive to low energy behaviour of the cross sections. 

8. CONCLUSION 

In a previous paper we have shown that new more accu- 
rate measurements of the CR antiproton flux pose a chal- 
lenge to existing CR propagation models. In particular, 
the antiproton flux and B/C ratio appear to be inconsis- 
tent with measurements when computed in standard dif- 
fusion/reaccelation models. In this paper we have demon- 
strated that this discrepancy can be resolved if some part 
of the CR that we measure near the earth consists of a 
"fresh" component accelerated in the LB. The indepen- 
dent evidence for SN activity in the solar vicinity in the 
last few Myr supports this idea. 

Combining the measurements of the antiproton flux and 
B/C ratio to fix the diffusion coefficient, we have been able 
to construct a model consistent with measurements of im- 
portant nuclei ratios in CR and derive elemental abun- 
dances in the LB. Calculated isotopic abundance distri- 
butions of Be and B are in perfect agreement with CR 
data. The abundances of three radioactive isotopes in CR, 
which are often used as "radioactive clocks" to determine 
the Galactic halo size, ""^"Be, ^^Al, ^^Cl are all consistent 
and indicate a halo size Zh ^ 4: kpc based on the most ac- 
curate data by ACE spacecraft. ^^Mn indicated a smaller 
halo, but this may be related to its half-life uncertainty 



and/or cross section errors. The derived fraction of LB 
component in CR is small compared to Galactic CR and 
has a steep spectrum with a cutoff above several hundred 
MeV/nucleon. Other experimental data (except maybe 
the overabundance of Sc, Ti, V) do not contradict this 
hypothesis. The derived source overabundance of sub-Fe 
elements in the LB may be caused by trivial uncertainties 
in the production cross sections, or could in principle arise 
from composition and/or evolution differences between the 
ISM in the LB and solar or Galactic average ISM. (LB may 
have evolved during the sun's lifetime of 4.5 Gyr.) This is 
suggested by the fact that the derived relative abundances 
of secondary elements in the LB sources are systematically 
larger than in the Galactic CR sources (see, e.g., derivation 
of interstellar ^"^C/^^C ratio, Savage et al. 2002). 

The production cross sections if measured accurately 
would help to distinguish between the different hypothe- 
ses; as of now, many important channels are not known 
accurately enough. Such cross section errors lead to errors 
in important isotopic ratios, which, in turn, are translated 
into errors in propagation parameters. In our treatment 
of Be and B production cross sections, as well as some 
isotopes of other elements, we use all available data and 
our own fits to them, which should be more accurate than 
semiempirical systematics by Webber et al. and Silberberg 
et al. 

The authors are grateful to M. Wiedenbeck for provid- 
ing the ACE isotopic abundances. I. V. M. is grateful 
to the Gamma Ray Group of the Max-Planck-Institut fiir 
extraterrestrische Physik, where a part of this work has 
been done, for hospitality. I. V. M. and S. G. M. acknowl- 
edge partial support from a NASA Astrophysics Theory 
Program grant. 



APPENDIX 

PRODUCTION CROSS SECTIONS OF ISOTOPES OF CARBON AND NITROGEN 

The production cross section of the most abundant ^^C isotope is one of the most poorly known. All the data available 
to us on its production by ^^O and nitrogen isotopes on protons are summarized in Table A3. The data include also 
production of ^^B, and ^^N, which decay to ^^C with branching ratio 0.98, and ^^O with branching ratio 0.12. 

Data on the production of ^^C are somewhat more extensive, but some important channels are not measured accurately 
enough. Most of the data available are summarized in Table A4. The data include also production of ^^N, which decays 
to ^^C with branching ratio 1. The production cross section of ^^B is very small, fractions of a mb. Fortunately there 
are data on production of ^"'N by protons on natural samples of oxygen and nitrogen (sec compilation by Mashnik et al. 
1998), which contain mostly ^^O and ^^N isotopes respectively. ^'^N production cross section data probably include also 
production of ^^O, but the latter cross section must be very small (0.17 mb at 2100 MeV/nucleon). 

The production cross section of the most abundant ^^N isotope is also poorly known. All the data available to us are 
summarized in Table A5. The compiled data include also production of ^^C, and ^^O, which decay to ^^N with branching 
ratio 1. 

The main contributor to the production cross section of ^'''N is ^^O (Table A6). The direct and indirect (via ^'"'O) 
production cross sections are almost equal. The channel ^^^O is well studied since there is a large amount of 

data obtained on natural sample of oxygen (see compilation by Mashnik et al. 1998). 

The cumulative (sum ov(^r all c;hanncls) production cross sections of carbon and nitrogen isotopes multiplied by the 
flux of the corresponding primary isotope in CR at 1 GeV/nucleon are shown in Fig. A15. The main contributor to the 
production of secondary carbon and nitrogen is ^^O, accounting for about 80% in case of nitrogen isotopes. However, in 
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the case of carbon, disintegration of ^^O gives only about 50% with an essential contribution from nitrogen isotopes (and 
i^C in case of ^^C). 

The contribution of '^'''N to the production of -'^''C is especially large (Table A4). This is based on only one experimental 
point which seems too large compared to the production cross sections on ^^N and ^^O. This indicates that the reason 
for the large fraction of ^^C in calculated CR abundances compared to the measurements (see Sections 2 and 6) may be 
errors in the cross sections. 

There are more examples of discrepancies in ^^C and ^'*N production cross sections (Tables A4, A5, shown in bold font). 
The production cross section of ^^C by ^^Ne at 580 MeV/nucleon differs significantly from that at 400 MeV/nucleon. The 
cross section of ^'^C by ^^Mg measured by the same group at 370 and 576 MeV/nucleon differs by a factor of 4 (6.3 mb vs. 
25 mb). A similar situation occurs with ^'^C production by ^^Ne at 400, 580, and 894 MeV/nucleon, and by ^^Mg at 371 
and 576 MeV/nucleon (3.5 mb vs. 9 mb). Fortunately, these latter cross sections do not contribute much to production 
of "'^^C and '^''N in CR, but these discrepancies indicate the degree of overall uncertainty in the production cross sections. 

Cross section errors in production of carbon may lead to errors in the B/C ratio, which, in turn, are translated into 
errors on the propagation parameters. Because the CR measurements are now rather accurate, the errors in the cross 
sections may cause many standard deviations when comparing the model calculations with CR data. 

It is clear that a more systematic approach to calculated cross sections is required such as using, e.g., evaluated cross 
sections in future similar work instead of only scarce experimental data or calculations by stand-alone nuclear reaction 
models or phenomenological systcmatics. Such evaluated data files (Mashnik ct al. 1998) have proved to be useful, 
e.g., to study the production of radioisotopes for medical and industrial applications using high power accelerators (Van 
Riper, Mashnik, & Wilson 2001). At present, neither available experimental data nor any of the current models or 
phenomenological systematics can be used alone to produce a reliable evaluated activation cross section library covering 
a wide range of target nuclides and incident energies. Instead, such an evaluated library may be created by constructing 
excitation functions using all available experimental data along with calculations employing some of the most reliable 
codes in the regions of targets and incident energies where they are most applicable. When there arc reliable experimental 
data, they, rather than model results, should be taken as the highest priority for the evaluation. The development of such 
evaluated data libraries for astrophysical applications is planned in the near future. 

Table A3 



Collection of C production cross section data on protons. 



Primary 


Secondary 


Energy, 


Cross section, 


Error, 




Nucleus 


Nucleus 


MeV/nucleon 


mbarn 


mbarn 


Reference 


14N 


120 


377 


56.90 


0.05*^ 


1 


14N 




516 


52.10 


0.05*^ 


2 


15N 




373 


30.00 


0.05^^ 


1 


160 




389 


33.90 


0.05*^ 


1 


160 




516 


33.60 


0.05*^ 


2 


160 




2100 


32.30 


4.80 


3 


160 


12b 


516 


1.10 


0.30^ 


2 


160 




2100 


1.45 


0.17 


3 




12n 


516 


1.10 


0.30^ 


2 


160 




516 


0.30 


0.30^ 


2 


160 




2100 


0.40 


0.07 


3 



* Relative error. 

''Relative error is shown as indicated in Webber et al. (1998a). 

References. — (1) Webber et al. (1998a); (2) Webber et al. (1990a); (3) Olson et al. (1983). 
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Table A4 



Collection of C production cross section data on protons. 



ST L LLLlciL y 


kJCCUllLldiy 


Ji/uergy, 




LLiL L DL , 




Nucleus 


Nucleus 


MeV/nucleon 


mbarn 


mbarn 


Reference 




13p 


oil 


7 fin 


u.uo 


1 
1 


IN 




010 


Q fin 


u.uo 


o 
z 


IN 




o to 


oD*oU 


u.uo 


1 
1 


16/-) 




ooy 


1 7 4n 




1 


i6n 




DID 


1 8 nn 


u.uo 


z 


16r\ 
\ I 






1 7 sn 

1 / .OU 


1 70 
± . / U 


Q 
O 


iMe 




^1^ 


1 ^ "^n 

lO. OU 


U.UO 


1 
1 


iNe 






1 ^ 7n 

10. ( u 


U.uo 


O 

z 


iNe 






1 7 sn 

1 ( .oU 


1 /in 


A 
41 


22 Mo 

iMe 




zini 


±O.OU 


u.uo 


1 
1 


22 Mp 




Ool 


91 Qn 


1 on 
1 .yu 


A 


22 Mo 

iMe 






1 Q nn 
ly.uu 


± .DU 




26 A/To- 
iVLg 




'?71 
Oil 


QaOU 


1 90 
1 .ZU 


/I 


iVLg 




^7R 




9 550 


/I 
41 




13N 




multi-data 




5 






377 


7.40 


o.w 


1 






516 


7.50 


0.10'^ 


2 


15N 




373 


4.30 


0.30'^ 


1 


nat0(160) 






multi-data 




5 


160 




389 


4.60 


0.20'^ 


1 


160 




516 


5.70 


0.20'' 


2 


160 




2100 


4.49 


0.46 


3 


20Ne 




414 


5.10 


0.20*^ 


1 


20Ne 




534 


4.10 


0.20^ 


2 


22Ne 




377 


0.50 


0.10 


4 


22Ne 




401 


1.60 


0.20*^ 


1 


22Ne 




581 


0.50 


0.10 


4 


22Ne 




894 


0.70 


0.20 


4 


24Mg 




610 


6.00 


0.20'^ 


2 


26 Mg 




371 


0.30 


0.10 


4 


26 Mg 




576 


0.10 


0.10 


4 



Note. — Discrepancy in the data is shown in bold. 
* Relative error. 

''Relative error is shown as indicated in Webber ct al. (1998a). 

References. — (1) Webber et al. (1998a); (2) Webber et al. (1990a); (3) Olson et al. (1983); (4) Chen et al. (1997); (5) a compilation by 
Mashnik et al. (1998). 
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Table A5 

Collection of production cross section data on protons. 



Primary 


Secondary 


Energy, 


Cross section, 


Error, 




Nucleus 


Nucleus 


MeV/nucleon 


mbarn 


mbarn 


Reference 


15N 




373 


27.60 


0.05'^ 


1 


160 


— 


389 


31.10 


0.05^^ 


1 


160 


— 


516 


31.00 


0.05'^ 


2 


160 


— 


2100 


31.00 


3.30 


3 


20Ne 




414 


25.80 


0.05^^ 


1 


22Ne 




401 


11.60 


O.IO'^ 


1 


15N 


140 


373 


10.30 


0.05^^ 


1 


160 




389 


1.70 


O.IO'^ 


1 


160 




516 


1.70 


O.IO'' 


2 


160 




2100 


3.69 


0.38 


3 


20Ne 




414 


2.20 


O.IO'' 


1 


20Ne 




534 


2.30 


0.10'^ 


2 


22Ne 




377 


8.10 


0.70 


4 


22Ne 




401 


7.70 


0.05^ 


1 


22Ne 




581 


10.20 


1.20 


4 


22Ne 




894 


8.60 


0.90 


4 


26Mg 




371 


3.50 


0.70 


4 


26Mg 




576 


9.00 


1.30 


4 


160 


140 


389 


1.30 


0.30^ 


1 


160 




516 


1.20 


0.30^^ 


2 


160 




2100 


0.75 


0.12 


3 


20Ne 




534 


1.00 


0.30'' 


2 


24Mg 




610 


1.50 


O.SO'' 


2 



Note. — Discrepancy in the data is shown in bold. 

^Relative error. 

''Relative error is shown as indicated in Webber et al. (1998a). 

References. — (1) Webber et al. (1998a); (2) Webber et al. (1990a); (3) Olson et al. (1983); (4) Chen et al. (1997). 
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Table A6 



Collection of °N production cross section data on protons. 



Primary 


Secondary 


Energy, 


Cross section, 


Error, 




Nucleus 


Nucleus 


MeV/nucleon 


mbarn 


mbarn 


Reference 


160 


15N 


389 


33.60 


0.05^ 


1 


160 


— 


516 


34.90 


0.05*^ 


3 


160 


— 


2100 


34.30 


3.30 


4 


20Ne 


— 


414 


24.00 


0.05=* 


1 


20Ne 




534 


27.80 


0.05'' 


3 


iNe 




o( ( 


on 







22Ne 




401 


32.90 


0.05'^ 


1 


22Ne 




581 


39.00 


2.50 


5 


22Ne 




894 


33.50 


2.10 


5 


24Mg 




610 


14.00 


O.IO*' 


3 


26 Mg 




371 


19.70 


2.10 


5 


26Mg 




576 


29.90 


3.00 


5 


22 Ne 


150 


377 


0.70 


0.20 


5 


22 Ne 





581 


0.80 


0.20 


5 


22 Ne 





894 


0.60 


0.10 


5 


26 Mg 


— 


371 


0.90 


0.30 


5 


26 Mg 


— 


576 


0.40 


0.30 


5 


nato(16o) 


150 




multi-data 




6 


160 




389 


30.70 


0.05^ 


1 


160 




516 


30.30 


0.05'" 


3 


160 




2100 


27.30 


2.60 


4 


20Ne 




414 


14.90 


0.10^ 


1 


20Ne 




534 


16.20 


0.10^ 


3 


21 Ne 




520 


7.80 


0.30'' 


2 


22Ne 




377 


2.00 


0.30 


5 


22Ne 





401 


1.60 


0.20'^ 


1 


22Ne 




581 


1.60 


0.30 


5 


22Ne 




894 


2.80 


0.40 


5 


22Na 




520 


10.90 


0.30*^ 


2 


23Na 




517 


11.10 


0.30^ 


2 


24Mg 




610 


8.60 


O.IO'' 


3 


25 Mg 




514 


6.00 


0.30*^ 


2 


26 Mg 




371 


1.10 


0.30 


5 


26Mg 




576 


0.20 


0.20 


5 


26 Al 




508 


4.10 


0.30^ 


2 


27 Al 




511 


7.10 


0.20'^ 


2 


28Si 




506 


5.60 


0.12^ 


2 


29Si 




508 


2.40 


0.30*^ 


2 



Note. — Discrepancy in the data is shown in bold. 



^Relative error. 

''Relative error is shown as indicated in Webber et al. (1998a). 

References. — (1) Webber ct al. (1998a); (2) Webber et al. (1998b); (3) Webber et al. (1990a); (4) Olson et al. (1983); (5) Chen et al. (1997); 
(6) a compilation by Mashnik et al. (1998). 
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Fig. A15. — Relative contributions of heavier isotopes to jModuction of carbon and nitrogen isotopes arc shown. These contributions are 
determined from cumulative cross sections of carbon and nitrogen isotopes weighted with the flux of corresponding primary isotope in CR at 
1 GeV/nucleon. The contributions of isotopes of a given element are indicated by separate lines. The scale gives the fraction of the C and N 
isotope produced by a given primary isotope. 
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